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ABSTRACT 


ANALYSIS AND PARTIAL DEVELOPMENT OF A SYSTEM 
FOR MEASUREMENT OF THE MOTIONS OF A SHIP 
MODEL IN SIX DEGREES OF FREEDOM 


by 
Richard 0, Haughey and Peter J. DeLaat, Jr. 


Submitted to the Department of Naval Architecture and Marine Engineering 
on 26 May 1958 in partial fulfillment of the requirements for the degree 
of Naval Engineer and the degree of Master of Science in Naval 
Architecture and Marine Engineering. 


The object of this investigation was to propose, analyse, and 
partially develop an angular and linear motion sensing system suitable 
for seaworthiness testing. It is anticipated that the principal use 
of the developed system will be in the measurement of model motions. 
In particular, ship model testing at the M.I.T. Towing Tank, and tanks 
of a similar character, was given full consideration, 


Emphasis was placed on the selection of transducers, a proper 
mounting scheme to isolate undesirable disturbances and cross coupling 
effects between the inputs, and a suitable method to resolve the 
Signals into independent components along the three chosen axes. 


The feasibility of a system employing the use of three linear 
accelerometers and three rate gyros is believed to be established. 
The use of such a system affords a continuous time base record of 
linear accelerations, velocities, and displacements, as well as a 
record of angular rates and motions. In view of the characteristics 
of such a sensing package, it is felt that weight and space limitations 
of model carried equipment are no longer a problem. 


The system, in its present state of development, will require 
additional study taking into consideration (a) refinement of the 
gravitational compensation system, and integrator balancing system, 
and (b) the laboratory mockup should include the use of rate gyros 
affording an accurate establishment of their performance 
characteristics. 


Thesis Supervisor: Y. T. Li 
Title: Associate Professor of Aeronautical Engineering 
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difference between 

yaw (rotation about the z axis) 

micro (10-6) 
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I. INTRODUCTION 


The problem to be treated in this investigation is the development 
of a system for accurately sensing the motions of a freely floated ship 
model. This definition implies one of the major constraints placed 
upon the system. That is to sense without inhibiting. The system will 
pive: full Sprttaatetice to the character of the M.I.T. and similar 
towing tanks, placing emphasis upon their size and facilities 
available. 

Recently there has been considerable interest, both analytical 
(43)" (25) and experimental (35) (39), in the problem of predicting 
the motions to be expected from a ship in rough seas, commonly 
termed seaworthiness evaluation. The most established means of 
accomplishing this goal is to make a study of model motions and then, 
by the use of Froude's Law of Geometric Similitude, to predict the 
desired ship motion. 

There have been numerous systems developed for the purpose of 
measuring either specific motions (7) (40) or combined motions of a 
model (39) (36). These systems have employed the use of motion 
pictures, synchro transmitters, gyro transmitters, mechanically 
linked potentiometer pickoffs, and electromagnetic coupling. Strobo- 
scopic (41) and flash-photography (42) methods have also been used. 
These systems suffer from one or more of the following disadvantages: 

a. specific motion sensing devices are not adaptable to the 
uses of the present problem. 


b. restraint of the model. 


* Refer to Literature Citations, Appendix E. 
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c. excessive size and weight. 

d. lack of a desirable level of precision. 

e, the necessity of time consuming data reduction, which has the 
twofold effect of preventing rapid data analysis, and preventing the 
elimination of secondary effects that destroy the precision of the 
test results. 

The system to be proposed uses a small sensing package, located 
at the center of gravity of the model, containing three linear 
accelerometers and three miniature rate gyros. Such a system will 
provide a continuous time base record of the three linear accelerations 
and their first and second integrals, as well as the three angular 
rates and their first integrals. There will be several major 
advantages to such a system. The overall weight of the sensing 
devices is a little less than a pound (this does not include a 
weight estimate of the mounting which they will require), and will 
occupy a volume approximately five inches by four inches, by four 
inches. The gyros mentioned in the earlier systems are of the free, 
or two gimbal, type and are not available in a miniaturized form. 
Consequently, their weight is of the order of six pounds (7), and the 
volume of a single unit is more than the entire package stated above. 
An inference may be drawn from the earlier portion of this paragraph 
that accelerations are considered more desirable than displacements. 
This is not consistent with the order of precedence established by 
the Thesis Advisor. However, it is felt that, by using the proposed 
system, it will be possible to obtain output signals of an accuracy 


somewhat superior to the presently available levels. (35) (39) 
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Ii. PROCEDURE 


A, Establishment of the System Specifications 


1. The system should be capable of measuring the displacements of 
the ship model in all six degrees of freedom. 

2 tae predicted range of accelerations and displacements to which 
the system will be subjected are tabulated on the following page in 
Table I. | 

3. The model carried components will not exceed five pounds 
maximum weight, and will not exceed the dimensions of four inches high, 
by four inches wide, by six to eight inches long. 

lL. The system will be made as accurate as possible, but the 
deviation should not in any case exceed 5 percent. 

5. The system should be capable of adhering to these limits of 
accuracy for a period from two to four minutes with some consideration 
being given to a ten minute maximum. 

6, It is intended that the overall system be capable of the 
measurement of both linear and angular quantities. In accordance with 
this EE cier, it was decided, after consultation with the Thesis 
Ravi cor, that the following precedence would be followed: 

a, Precedence of signal resolution 
displacement 
velocity 
acceleration 

b. Precedence of motion measurement 
roll 
pitch 
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surge 


sidle 


yaw 


B. Factors Considered in the Selection of Relative Motion Receivers (1) 


1. 


Is 


Type of Input. 

linear 

angular 

Input range. 

Maximum uncertainty in the low input magnitude range. 
Maximum uncertainty in the high input magnitude range. 
Dynamic characteristics. 

Availability. 

Size and shape of space occupied. 

Weight. 

Useful life and reliability. 

Maintenance. 


Cost. 


C. Factors Considered in Relation to the Over-all System Performance (1) 


ibe 


Sensitivity to electrical interference. 


2. overvices required for operation. 
3. Ease of operation. 
4, Size and weight of complete system. 
d¢ Reliability. 
6. Cost. 
D. Possible es of Sensing Elements Considered 
1. Linear accelerometers. 
2. Angular accelerometers. 
3, Rate gyros. 











4, Integrating gyros. 


5. Free gyros. 


E. Possible Combinations of Sensing Elements to Measure Displacements 
in Six Degrees of Freedom 


1. Six linear accelerometers. 

2. Thrée linear accelerometers and three integrating gryos. 

3. Three linear accelerometers and three rate gyros. 

4. Three linear accelerometers and two free gyros. 

5. One linear accelerometer capable of sensing the three 
mutually perpendicular linear accelerations, and two free 
gyros. 

6. Three linear accelerometers and three angular accelerometers. 

7. Six gyros with adjoining masses on beams. 

8. Three linear accelerometers capable of simultaneously sensing 
linear accelerations along two mutually perpendicular axes. | 


F, A General Discussion of Some Factors Relating to the System Design 





and the Selection of Sensing Elements 


The proposed sensing system is to be mounted on an unstable plat-= 
form, which will be simultaneously rotating and translating, This mode 
of excitation is to the best of the authors’ knowledge a new application 
for the proposed linear accelerometers, 

It was established, after consultation with the Thesis Advisor, 
that the system should sense motions with respect to a coordinate 
system which is fixed with respect to the ship model, as opposed to 
an inertial coordinate system. This is a rather obvious conclusion, 
since it is only those motions along ship axes which are consistent with 
the stated objectives of the thesis. The thought may occur Hove, 


however, that it would have been considerably simpler to consider an 


== 





inertial coordinate system, and to make use of motion sensing systems 
of a similar character already known to exist in the pilotless aircraft, 
missile, or rocket fields. This is an accurate, but somewhat specious 
thought. The existing systems are quite bulky (22), in relation to the 
stated specifications of the desired ss and would require a 
coordinate transformation. This latter requirement would have several 
disadvantages. (a) It would be complex. (b) If this coordinate trans-~ 
formation were feasible electronically it would probably be both 
expensive and, within the details of its own design, quite complex. 
(c) If it were not feasible electronically it would require an amount 
of data reduction, which would be 4 serious disadvantage in relation 
to existing model sensing systems. In this same vein it may be 
mentioned that the accelerometer package of an inertial guidance 
system, whether of the pendulous or single axis linear accelerometer 
type, is unsuited to the use of this design. There are a number of 
reasons for this. Some of the more obvious are (a) it is mounted on 
what amounts to a stable platform (23), and consequently suffers from 
the coordinate transformation difficulties mentioned above, and (b) in 
the process of arriving at the stable platform it makes use of 
integrating gyros, which are themselves undesirable pieces of equipment, 
as will be demonstrated below. 

In accordance with the intentions of the Preliminary Report dated 
22 January 1958 which states "—--the first portion of the thesis will 
be devoted to the study, and analysis of various systems. Upon com= 
pletion of this portion, a single system will be selected for further 
development. If the components of the selected system are available 
the thesis will proceed with the construction, calibration, and testing 


of this system. If the components are not available an evaluation of 
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the possibility of the construction of the desired components will be 
made,=--="" section "E" of this procedure was prepared with the specific 
intention of not limiting consideration of the possible combinations of 
transducers to known and available commercial models. An effort was 
made rather to consider the physical characteristics of the desired 
sensing system in light of the stated specifications, and while 
generally channeling the creative thought to those first principles 
utilized by successful transducing devices, an attempt was made to 
propose an instrument which might be better suited to the particular 
design problem under consideration. Reference is being made here, in 
particular, to linear accelerometers. 

It may be stated at this point that the same device which is used 
for the measurement of linear accelerations may also be used for the 
direct measurement of displacements (a seismograph), provided the 
natural frequency of the instrument is properly chosen. In light of 
this fact, the question might be very reasonably asked, if it has not 
occurred already, "Why were not displacement meters, rather than 
linear accelerometers, chosen as transducers in accordance with the 
precedence of signal resolution specified on page 37?" The answer 
lies in the fact that it is considerably easier (though by no means 
a small task) to double integrate than it is to accurately carry out 
a process of double differentiation (5). A somewhat contrasting, but 
nevertheless consistent opinion concerning the choice of linear 
accelerometers is contained in the following quotation; (23) 

"A satisfactory electrical integrator or differentiator within the 
optimum frequency of a vibrometer or accelerometer is relatively easy 
to obtain. For this reason an evaluation of the relative merit of an 


accelerometer system against a vibrometer system can be made, based 
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ae 


upon the property of the pickup alone. Generally speaking, an 
accelerometer is more desirable than a vibrometer. Firstly, for low 
frequency application of less than 5 cps, there are few satisfactory 
vibrometers available. Anything below 1 cps is almost nonexistent. 
----=- The construction of an accelerometer tends to be much smaller 
than a vibrofeter. This is especially true if a large amplitude of 
vibration is to be measured. Because no suspension problem is in- 
aie an accelerometer tends to be simpler and should therefore be 
less expensive. For this reason, the vibrometer as an instrument has 
little future, although there still seems to be a market for it." 

One more question of a general nature which may be asked at this 
time is, "Doesn't the vessel in fact roll about the virtual center of 
gravity (24) instead of the actual center of gravity, and since the 
linear accelerometers must be placed at this center of gravity in 
order to reduce the somewhat considerable errors (developed below) 
due to angular motions, the system might well prove to be impractical?" 
Here again it is rather easy to confuse the elements of this design 
problem with those of existing analyses. The concept of the virtual 
center of gravity results from a consideration of the virtual mass of 
a body used for the determination of the inertia forces experienced 
by a model in a perfect fluid. These inertia forces are then used 
for the analytical prediction of the motions of a ship. (25) The 
concept of added mass together with that of the virtual center of 
gravity is then seen to be an analytical concept employed for the 
successful prediction of ship motions. The proposed sensing system 
is not concerned with the prediction 2 but the measurement of the 
model motions, and consequently need not concern itself with the 


concept of the virtual center of gravity or its movements. 
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G. Specific Analysis of the Possible es of Sensing Elements Considered 
The analysis contained below is somewhat abridged from its original 
form. It is felt that it would be redundant to repeat in the case of 
each instrument the numerous factors mentioned under paragraphs "B" and 
"Cc" above. Rather, reference will be made to the particular manufac- 
turer's catatogue from which the desired information may be obtained. 
In other words, the intention is to present those elements of the 
analysis which are more closely related to the particular design under 
consideration and, by so doing, to prevent occlusion of some of its 
more salient features by an unnecessary amount of detail. 
1. Linear Accelerometers: 

(a) Linear accelerometers, regardless of type, are essentially 
simple, seismic, second order, single-degree-of-freedom systems. 
Table II (19) lists some of the more prevalent types and presents 
their general characteristics. 

(b>) The signal from the linear accelerometers will require 
double integration before displacements are obtained. 

(c) It was originally proposed to use sum and difference 
operations on the output of the two linear accelerometers in the 
manner presented in Figure I, with the instruments displaced from 
one another along a ship's axis. Figure II” demonstrates that at 
the maximum displacement allowed by the proposed system specifications 
(two inches either side of the system center of gravity), the tangential 
acceleration due to the largest predicted angular acceleration (roll) 
is only .96 fps“. The order of magnitude of this signal may be 


contrasted visually with those of the predicted linear signals (surge, 





* Por sample calculations see Appendix D. 
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heave, and sidle). It was felt that such a low level signal might well 
be obscured by errors due to non-linearity alone, and in view of the 
motion measurement precedence, previously established, the use of six 
linear accelerometers as possible angular sensing devices was effective- 
ly ruled out. 

(a) At this point in the procedure it becomes clear that the 
sensing of angular accelerations by the linear accelerometers, either 
normal or tangential, is an undesirable occurrence, An ideal design 
for the purpose of reducing angular acceleration effects would be one 
linear accelerometer capable of sensing the three mutually perpendi- 
cular linear accelerations, to be placed at the center of gravity of 
the model. Consistent with the general intentions of the er some 
thought was given to the design of such a device. This avenue of in= 
vestigation was concluded with the observation that since existing 
accelerometers cannot be designed to within a closer tolerance than 
0,02 g per g for response to transverse accelerations, a continued 
pursuit of this study would probably produce an instrument whose 
characteristic in this respect is somewhat in excess of the above, 
and consequently unsuitable for the purposes of the proposed system. 

2. Angular Accelerometers: 

Angular accelerometers possess response specifications (11) 
suited to the requirements of the proposed system. However, they are 
unsatisfactory from the point of view of weight and space, and con= 
sequently were not selected for further consideration. 

3. Rate Gyros, Integrating Gyros, and Free Gyros: 

In contrast to the individual treatment of the linear and 

angular accelerometers, it is felt that there is little to be gained 


from such a breakdown of the analysis of the various types of gyros. 
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While there are a number of considerations to be applied in the selection 
of the type of gyro to be used for the measurement of a particular input, 
the choice was limited in this case by the size specifications of the 
transducer package. The result was the selection of the miniature rate 
gyro. (8), (9), (10). An evaluation of the errors anticipated from the 
use of this instrument is included in the Discussion of Results. 
H. Analysis of Linear Accelerometer Errors Due to Rotational Motions of 
the Model 

In view of the character of the acceleration sensing devices, it 
is apparent that any displacement of their sensitive masses from the 
center of gravity of the model will result in an output signal due to 
the algebraic sum of those components of normal and tangential 
accelerations coincident with their sensitive axes. At this point in 
the procedure it was realized that the most direct optimization process 
for minimizing the errors resulting from this effect would be to 
develop three sets of curves per axis for various positions of an 
accelerometer along any given axes. These curves are plotted on 
Figure XII. Sample calculations for these curves are shown in 


Appendix D. 


I, Choice of the Specific Linear Accelerometer Models and Their 
Proposed Positional Arrangement 


In response to correspondence with the Statham Laboratories, 
specifications (3) and data (12) were obtained for the purpose of 
selecting particular accelerometer models. Heavy paper mockups of the 
various models were constructed, facilitating a realistic determination 
of the accelerometer package size. Numerous configurations of the 
models were tried in order to arrive at an optimum arrangement. This 


procedure considered, in addition to those factors mentioned above, 


15s 





(a) Figure XII, (b) the precedence of motion measurement, and (c) a 
desirable damping coefficient. In view of the system specifications 
and (12) there is good reason to believe that a damping coefficient 

of something less than .7 would be desirable. However, this nominal 
value was adhered to with the intention of not limiting the usefulness, 
or number of possible applications, of the proposed instruments. The 
selected instruments and their configuration are shown in Figure XIII. 
J. Choice of the Specific Rate Gyro 

The selection of satietactory rate gyro models has been mentioned 
on page 15. This paragraph is entered primarily for the purpose of 
continuity and to mention that it was determined, after their 
selection, that the rate gyros could not be delivered within a period 
of time that would be satisfactory for the Thesis. 

K, Static Calibration of the Linear Accelerometers 

The accelerometers were simultaneously calibrated and examined 
for linearity and i oteresi effects by the simple procedure of rotating 
them, at fixed intervals, through 360 degrees in the earth's gravi- 
tational field. Tabulated data for these tests is contained in 
Appendix C. The results are shown on Figures XIV through XIX. 

The calibrating setup was composed of equipment listed in 
Appendix E, Prior to commencement of calibration, the longitudinal 
axis of the Hardinge dividing head was aligned horizontally by the use 
of precision levels, Then the accelerometer models were placed in a 
4, jaw chuck on the dividing head and their cases aligned with the long 
dimension of the instrument vertical, and the short dimension brought 
into the horizontal plane by means of adjusting and zero setting the 
calibrated face of the dividing head. With the instruments thus 


aligned, they were provided with battery excitation and their output 


SiGe 





recorded at the various angles of inclination by means of a precision 
eee Poet stat in conjunction with an auxiliary, or external galvan- 
ometer. (38) 
L. Consideration of the Instruments Response to Transverse Accelerations 

The manufacturer's specifications (3) for these instruments states 
that they have a response to the vector sum of those accelerations in a 
plane cr their sensitive axis of .,02 g per g. Figure XX was 
prepared with the intention of evaluating the predicted variation of 
this error and is believed Pe be self explanatory. 
M. Gravitational Error Analysis 

This analysis was somewhat lengthy and consisted of graphical 
presentations of a number of somewhat complex expressions for the 
purpose of arriving at a simple but satisfactory expression which 
might lend itself to electronic reproduction. The procedure followed 
here was to (a) develop curves defining the actual variation of the 
gravitational error present during the various motion measurements, 
(b) consider approximation of the error by making use of a Taylor's 
series expansion, and (c) consider approximation of the error by the 
use of trigonometric functions. Results of this procedure are shown 
on Figures XXI through XXIV. Derivation of equations used is contained 
in Appendix B, and a summary of calculations in Appendix C. 
N. The Proposed Sensing System 

The proposed system, Figure XXVIII, was based on all of the fore- 
going procedure in addition to the following considerations: 

(a) signal attenuation due to integration.” 


(b) trigonometric function generation by means of the 


* See Appendix D. 
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"Sine-Cosine Operator." (26) 

(c) trigonometric function generation by means of servo driven 
sine-cosine potentiometers. (21), (27), (31), and (32) 

(d) characteristics of, and type of recording system. (33) 

O. Laboratory Procedure 

After completing the proposal of the desired sensing system, an 
effort was cal to ascertain its feasibility by submitting a number 
of its components to laboratory test. 

Prior to the commencement of the tests, however, it was necessary 
to assemble the Double Integrating Summing Amplifier shown pictorially 
in Figure III, along with a suitable calibrating platform shown pic- 
torially in Figure IV. The Double Integrating Summing Amplifier made 
use of two Model USA-3, Philbrick operational amplifiers. (20) A 
schematic presentation of a single stage is shown in Figure V, and a 
derivation of its theoretical frequency response is included in 
Appendix B. The calibrating table is a relatively simple device which 
allows the accelerometer to be subjected, simultaneously, to a static 
gravitational acceleration, at various angles of inclination, while 
experiencing dynamic accelerations consistent with the predicted am= 
plitudes stated in the system specifications. 

The tests which were carried out are as follows: 

1, Calibration of the Sine-Cosine Potentiometer. 

The procedure employed for the static calibration of the sine- 
cosine potentiometer was quite similar to that employed for the static 
calibration of the linear accelerometers. The equipment used is listed 
in Appendix E, a sample calculation for the voltage divider circuit is 
included in Appendix D, and a schematic diagram is included in Figure VI. 


The shaft of the potentiometer was secured in the jaws of the dividing 
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head while taking care to support the potentiometer housing in such a 
manner as to prevent wobbling of the shaft while rotating it through 
360 degrees. 

2. Ability of the Sine-Cosine Potentiometer to Null Gravitational 
Brrors. 

With the accelerometer and potentiometer mounted as shown in 
Figure IV, the output signals of the two instruments were recorded in 
Table XI , and compared graphically in Figure XXX, 

The essential steps carried out during this test are (a) the 
outputs of the two instruments were established at zero with the table 
horizontal, (b) the table was rotated to an angle of forty-five degrees, 
and the accelerometer output was measured, (c) by varying the decade 
resistor, shown in the schematic diagram Figure VII, the magnitude of 
the voltage divider pickoff was made the same as that of the 
accelerometer output, and (d) with the system thus adjusted the table 
was rotated back into the horizontal, in five degree steps, while 
recording the output of the two instruments. 

3. Double Integrator Response to a Constant Displacement Input. 

Prior to the commencement of this test, the output of the 
accelerometer was established at zero with the table horizontal. A 
schematic diagram of the circuit components and their connection is 
shown in Figure VIII. A picture of the components is shown in Figure 
IX, and an equipment list is included in Appendix E, 

The full travel of the Scotch Yoke, Figure X, was measured with a 
steel rule and found to be 4.763 inches. The initial speed of the yoke 
was determined to be 76 rpm by use of a stop watch. 

With the equipment set up as described aoueee the rpm of the 


Scotch Yoke was varied over the frequency range shown in Figure XXXI. 
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The different amplitudes thus obtained were normalized with respect 
to that for 76 rpm and plotted on Figure XXXII. 

The frequency response, Figure XXXTIT, of the double integrator 
was calculated from this same recording by scaling the output in 
millivolts from the Sanborn recorder data. The acceleration signal 
was calculated for the known acceleration amplitudes of the 
Brees ceiicter be usiie the calibration factor for that instrument. 

The system gain was then divided by 100 (the Video amplifier gain) (28) 
in order to obtain the gain of the double integrator. 

1, Establishment of the Drift Characteristic of the Double 
Integrating Summing Amplifier. 

The experimental setup for this test was the same as that des- 
cribed for test number three, There is little need here for amplifying 
detail except to point out that the test was run for ten minutes, and 
the combination of yoke speed (76 rpm) and displacement (4.763 inches) 
results in an acceleration of 12.5 fps< which is the predicted maximum 
acceleration. 

5. Dynamic Gravitational Error Compensation of the Linear 
Accelerometer Signal at Fixed Angles of Inclination. 

The experimental setup was once again quite similar to that 
described for test number three, only this time it included the 
Sine-cosine potentiometer circuitry, which may be seen by an in- 
spection of Figure XI. 

The essential steps carried out during this test are (a) the 
outputs of the accelerometer, and the sine-cosine potentiometer were 
established at zero with the table horizontal, (b) the table was 
rotated to an angle of forty-five degrees, and the accelerometer output 


was measured, (c) by varying the decade resistor, the magnitude of the 


28. 
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voltage divider pickoff was made the same as that of the accelerometer 
output, and (d) with the system thus adjusted the Scotch Yoke was 
activated, subjecting the accelerometer to a known acceleration 
amplitude of 12.5 fps“. It was intended to rotate the calibrating 
table back into the horizontal plane, at five degree intervals, while 
noting the system output. However, the test results were unsuccessful 


and, due to lack of time, this testing procedure was discontinued. 
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Iii. RESULTS 


Sensing Elements 


a. Linear Motions: Linear Accelerometer, unbonded strain gage type. 
Statham Laboratories, Inc. 
Natural Freq.: 21-40 cps. 


Transduction: Resistive, complete balanced bridge 
% 


Excitation: D.C. or A.C. 9-11 volts 
Output : 17-35 mv per g. 
Damping: 0.7 (+ 0.1) of critical at room temp. 


Ambient temp. limits: -40° to + 150°F. 

Max. allowable static acceleration: 3 times rated range 
Response to transverse acceleration: .02 g per g. 
Non-linearity and hysteresis: 1% of full scale excursion. 
Weight: 35 grams to 13.5 ounces 


Size: 33" x 2x yin or smaller. 


Soy 





b. Angular Motions: 


Full scale rate: 
Angular momentum: 
Undamped natural freq: 


Damping Ratio: 


Excitation: 


Full scale gimbal 
deflection: 


Acceleration time: 
Full scale output: 
Sensitivity: 

Threshold & Resolution: 
Inertia Ratio: 


Type of Pickoff: 
Linearity: 


Ambient Temp. Range: 


Shock or Linear 
Acceleration: 


Linear Vibration: 


Size: 


Weight: 


Miniature Rate Gyroscopes. 


Manufacturer's Data 


Minneapolis Honeywell 
GN-40 or 


L00°/sec. 

38,600 gm-cm@/sec. 
76 cps. 

O%6 


6.3 v, 400 cps, 2 phase 
26 v, 400 cps, 2 phase 


Bao 


SaSee, 

5.6 v. (RMS) 

14 mv/geg/sec. 

Ol degrees/sec. 
~0O00L5 sec. 
Variable Reluctance 


0.1% to 4 range 
2% to full range 


~67° to + 185°F 
Up to 100 g. 


Up to 10 g,O to 2000 cps. 


1" diam. x 2£" long 


3.8 ounces 


aaoe 


Ui e. lime Corp. 
UST #400 


+ 1,00°/sec. 
32,600 gm-cm@/sec. 
88 cps. 

Of 7 onretorsurt 


6.3 v, 400 cps, 2 phase 
26 v, 400 cps, 2 phase 


20g 


15 sec. 
5.6 v. (RMS) 


14 mv/dep/sec. 

.02 degrees/sec. 
e000 sec. 

Variable Reluctance 


0.1% to 4 range 
1% to full range 


-67° to 212F. 


>200 g. 


25-30 g. 


94" diam. x 2 3/16" 
long. 


3.5 ounces. 
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Due to Accelerometer Displacement From the Model C.G. 


FIGURE XIl 


Linear Accelerometer Direct Axis Error Due 
To Anqular Velocities and Accelerations 





(For Tabulation of Data. See Appendix B) 


NOTATION: 
Accelerometer 


Axis along which the 
Accerometer C.G. is 
Displaced from the 





Distance in inches 


E336 





FIGURE Xi U1 


Accelerometer Arrangement 


Ref.: STATHAM INSTALLATION DWGS 6620 and 6623 


Note: All Dimensions are in Inches 


TOP VIEW Full Scale 
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FIGURE XX 


Maximum Error in Statham Linear Accelerometers 
Due to Linear Accelerations Normal to 
The Indicated Sensitive. Axis* 
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Transverse Acceleration in FPS* 
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Any acceleration normal to the indicated axis 


Vector sum of A (heave or surge) and. maximum sidle signal 
_ Va2 2 
B = VA“ + 36 (ay wax ) 


Vector sum of A (surge or sidle) and maximum heave signal 
c =Va2 + 100 


Vector sum of A (heave or sidle) and maximum surge signal 


6 FPS 


10 FPS2) 


(2 -max 


D =\/A2 + 156.2 (apnay = 12.5 FPS) 


Ref. (3) 
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GRAVITATIONAL ERROR 
IN HEAVE FOR VARIOUS 
ANGLES OF ROLL AND PITCH 
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Gravitational ‘Error in Surge for 
Various Angles of Pitch and Yaw 
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Calibration of Sine-Cosine Potentiometer 








FIGURE XXX 


Ability of the Sine-Cosine Potentiometer 


To Null Gravitational Errors 
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Accelerometer Inclination 
(Degrees) 


‘Accelerometer: Statham A6a=1.5-350; Ser. No. 3023 
Excitation: 8.4 Volts (0.C.) 


Sine-Cosine Potentiometer: Gamewel | RL=-14MS: Ser. No. 674 
Excitation: 7.8 volts (D.C.) 
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IV. DISCUSSION OF RESULTS 


A, Selection of Sensing Elements 


The sensing elements chosen fall well within the space and weight 
restrictions of the model carried equipment. It is felt that these 
instruments represent the optimum selection of commercially available 
components. The overall system is not the least complex that could 
have been designed, but the space and weight restrictions of the 
model carried components have imposed this compromise. 

The time required for delivery did not permit the purchase of 
the rate gyros. However, a brief analysis of possible gyro errors 
shows the following: 

1. Hysteresis is .05% of full scale or 0.2° per second. The 
smallest design magnitude of angular velocity is that for 
yaw (39.2° per second). Therefore, the largest possible 
error is only .51%. 

2. Linearity (to 200° per second) is 0.1% of full scale or 
0.4° per second. The largest possible error (yaw) would 
be 1.02%. 

3. Variation of sensitivity with temperature is .02% per °c. 

A 50° variation would introduce only 1% error and temper- 
ature variations of this nature are not anticipated. 

In view of the above, it is felt that either the Golden Gnat 40 

(10) or the UST #400 (9) would provide signals well within the error 
limitations. Significant errors, if present, would arise from other 
components in the particular system, 

U. S. Time Corporation has advertised telemetering packages with 


D.C. output and 3 axis control packages (9). It is recommended that 
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these be investigated further and perhaps applied to this system. 

From an inspection of Table II, it is apparent that the only 
linear accelerometers which meet stated system specifications are the 
strain guage type and the differential transformer type. The unbonded 
strain guage type (3) was selected in lieu of the differential trans- 
former after a general consideration of the type of power supply 
required ey ithe two instruments. It was felt that the 400 cps 
carrier supply required for the differential transformer would be 
subject to the usual frailties of regulated power supplies (17), and 
might unnecessarily complicate the overall system. 

Figure XII played a very important role in choosing specific 
accelerometers as well as locating them within the package. This 
method of analysis is rather unique and useful. Having once es- 
tablished the anticipated angular accelerations, a simple calculation 
(page 86) will define this plot and show at a glance the best 
arrangement along the axes for accelerometers in the package, the 
minimum errors being governed only by the physical dimensions of the 
accelerometers. Figure XIII shows the result of just such an analysis. 
Other combinations of accelerometers approached this minimum error but 
were discarded because the output range of the sensing elements was 
smaller, The A5a~1.5~350 accelerometer senses heave, the priority 
linear signal. It is therefore located at the center of gravity to 
eliminate all errors due to roll and pitch accelerations. It is also 
the best choice in shape to allow the other two to be located as close 
as possible to the center of gravity of the model. 

As the model rotates, the accelerometers sense that component of 
the gravitational acceleration parallel to their sensitive axes. 


Since the heave accelerometer senses 1 gravity when in its neutral 
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position, it must be capable of measuring this signal as well as heave. 
A 1.5 g instrument was chosen. Instruments are available for vertical 
measurement that use an internal spring to bias them to the neutral 
position. This would be of no benefit in this system, because the 
instrument would simply measure gravity in the opposite direction as 
the model rolled. Likewise, the sidle accelerometer, also affected 

by roll, was fence to be 1.5 g. To allow some degree of standard- 
ization of components, the same sensing magnitude was chosen for the 
surge accelerometer. 

It is felt that perhaps the specification of a maximum roll 
angle of 60° is excessive. If 45° were specified, the maximum 
variation in the gravity signal would be about 0.7 g. Since the 
maximum heave signal specified is approximately .31 g, al g in- 
strument could be used with very infrequent possibility of exceeding 
the rated accelerometer range. This would provide two important 
advantages over the selected instruments: 

(a) the output sensitivity would be increased by about fifty 

percent, and 

(b) the non-linearity and hysteresis error would be reduced 

by approximately the same percentage, since it is rated 
at 1% of full scale output and we are now reducing full 
scale output as well as restricting the signal to something 
less than before. 
This improvement would be noticed most in the sidle accelerometer, 
for its desired signal is specified as only half that of the other 
two accelerometers. 


B, Calibration and Error Analysis of the Linear Accelerometers 


The calibration factors, as determined experimentally, are in 
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very close agreement with those furnished by the manufacturer. It 
should be noted that the calibration curves do not pass through the 
origin on Figures XVI and XVIII. This is due to the hysteresis of 

the instruments, and the discrepancy could probably be reduced by 
lowering the damping factors of these instruments and making it easier 


to zero the instruments when installing them in the model. However, 


¢ 


if the instruments are installed as they now exist, proper zeroing 
procedures should reduce this discrepancy to negligible proportions. 

The non-linearity and hysteresis errors are well within the 
manufacturer's specifications and should cause no trouble in data 
recording. On the surface, the errors appear to be greater than 
allowed by the system specifications, especially in the surge and 
sidle accelerometers. However, it must be remembered that in the 
actual model, a full excursion will not be experienced. The surge 
accelerometer will experience gravity due to pitching, but the maxi- 
mum pitch is only 10°. This means the maximum gravity signal is only 
about one-half of gravity, thus restricting the non-linearity and 
hysteresis error to approximately one percent of full scale surge. 
The sidle accelerometer will experience about one gravity (due to 
roll and sidle) at its maximum signal. The error introduced here will 
be about 63% of the maximum sidle signal which means an accuracy 
within one-tenth of an inch for sidle displacement. It is anticipated 
that sixty degree rolls will be very rare and therefore, the error 
much reduced. 

It is recommended that the excitation voltages be as high as 
possible (up to 9 volts maximum) in order to insure maximum accuracy 
and resolution (23). A large capacity battery source should be used 


to insure a constant level of voltage during model runs, and to 
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eliminate the presence of ripple ordinarily found in power supplies. 

Figure XX is a plot showing the range of errors in acceleration 
due to linear accelerations normal to the sensitive axis. This range 
is defined as the area between curve A and any one of the other three 
curves and extends out to the line of maximm signal consistent with 
the signal assumed as A. For example, the maximum error in surge due 
to heave and sidle is .234 fps* or 1.87%. This chart provides a rapid 
means of determining the maximum error along any axis due to acceler- 
ations along the other two. The chart is based on the maximum 
sensitivity to transverse accelerations as specified by the manu- 
facturer. It has been predicated by the thesis supervisor that this 
error is caused by misalignment between the sensing element and its 
case and can be eliminated (or reduced to a negligible quantity) by 
aligning the mass axis, rather than the case axis, with the model. 
This is done by rotating the accelerometer until its signal is nulled 
for all attitudes of the accelerometer when rotated about its sensitive 
axis. The accuracy of such a procedure, however, is limited by the 
hysteresis of the instrument. Since this particular error is an un- 
certainty error, it will be disregardédin the consideration of com- 
pensation systems for error signals. 

The major error signal is that due to gravity. For example, the 
heave accelerometer senses a full gravity when vertical and may sense 
as little as a half g when rolled to 60 degrees. This is almost twice 
the signal it is desired to measure and therefore must be eliminated 
by a compensation system. Exact expression for the gravity error due 
to rotational motions have been developed in Appendix A for each of 


the three linear accelerometers. 
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(Heave) e=g pu cos @ copy 


ay 
(cos“d + cos’ ¥ - cos“ cos“ )F wy) 
(Surge) e=g sin ¥cos 6 as 
(cos*O + sin? Y - sin® Y cos*6) 
(Sidle) e=g sin § cos 6 - 
(cos“O + sin2g - sin~ g cos 6)2 3 


The sensing system generates the three angles, so all that is left 
for the compensation system to do is to resolve these angles into 
usable form. Results of analysis of the heave error are shown in 
detail and a summary comparison is made for the other two linear motions. 

Figures XXI and XXII show that the roll compensation system must 
include inputs from both roll and pitch to keep the error within 
acceptable limits. 

A Taylor series approximation of the expression (Figure XXTIT) 
showed that a fourth order equation would be required to keep the error 
within acceptable limits. Therefore, this method of analysis was 
abandoned in favor of a more direct type of trigonometric function 
generator. 

Figure XXIV shows the results of attempts to simplify equation (1) 
by fairly simple approximations. The expression"e = g cos J cos ¥" 
shows great promise, since the error at simultaneous extreme roll and 
pitch is somewhat less than two percent. 

Taylor series expansions for the surge and sidle expressions also 
required an equation of too high an order to be feasible (Appendix A). 

Simplification in terms of trigonometric functions provided ex- 


pressions somewhat easier to generate and the errors from these 
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approximations are considered negligible (Figures XXV and XXVI). 
C. Proposed Gravity Error Compensation System 

The proposed gravity compensation system is shown in Figure XXVII. 
Two systems were deemed possible for generating the trigonometric 
functions: 

(a) servo driven sine-cosine potentiometers (21), (27) 

(b) elecfronic sine-cosine operator (26). 
Servo drives are used in most low frequency analogue computer work, 
since their response to low frequency signals is quite satisfactory. 
The sine-cosine operator is intended for use in high frequency 
repetitive systems. The high cost of such electronic units, along 
with the fact that either the sine or the cosine, but not both, is 
generated per unit, led to the selection of the servo system for 
generation of the gravity compensation signals. 

No experimental work has been done with the proposed system, but 
a Gamewell (21) potentiometer was used in the laboratory simulation 
model (Figure XI) and some difficulty was experienced in nulling the 
gravitational signal. The sine-cosine operator was therefore mentioned 
above to call attention to it as a possible substitute should further 
experimentation prove the potentiometer unsatisfactory. 
D. Proposed System 

The complete system, exclusive of power supplies, is shown in 
block diagram form in Figure XXVIII. A mockup and partial analysis 
of the linear displacement system for one axis has been made and 
results (to be discussed later) look very promising. 

The only equipment to be carried in the model is the package unit 
consisting of three linear accelerometers and three rate gyros. These 


components are to be placed in a small box known as the sensing unit. 
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The center of gravity and axes of the unit will be clearly marked on 
the unit case and should be placed exactly at the model center of 
gravity and the axes aligned with the model axes with extreme care. 
Failure to comply with this precaution will introduce serious errors 
into the linear acceleration sensing system. The accelerometers will 
be mounted as shown in Figure XIII. 

The D.C, amplifiers (28) in the linear accelerometer circuits 
were originally included to bring the accelerometer outputs to satis-—- 
factory levels for use in the inves ra vim amplifiers. However, they 
may prove to be unnecessary since the double integrating amplifier 
as it now exists (Figure V) has a gain built into it. 

The double integrating amplifier is based on Philbrick USA-3 
operational amplifiers (20)", It works very well, once properly 
balanced, but balancing is a tedious process, and further refinement 
is necessary before the unit can be used effectively at the model basin. 
A basic refinement would be the replacement of the present balance 
potentiometer by a ten turn potentiometer with plus and minus voltages 
introduced from the two ends rather than with the switch presently 
installed. The magnitude of the balance voltage available should be 
increased to allow for balancing the amplifier over a greater range 
of unbalance, 


E. Laboratory Data 


The calibration curve for the sine-cosine potentiometer (21) is 


Supposedly, there are several good chopper stabilized D.C. operational 
amplifiers on the market, but no satisfactory ones are available at 
M.I.T. The Philbrick K3-J integrating amplifier (29), REAC integrating 
amplifier (30), and the DACL computer integrating amplifiers (no 
publications in existence) all have noise levels referred to the input 
that are in excess of those acceptable for this system. 
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quite erratic. Errors shown are a percent of the sine (or of gravity) 
and therefore would be higher when referred to specific model motions. 
The potentiometer used in the laboratory mockup is a rather old one 
and may therefore be worn. The fact that practically all errors are 
in the same direction may raise some doubt as to their validity, but 
in fact this is of little consequence, for it merely indicates poor 
zeroing. The erratic nature of the error is of most importance. The 
method of calibration is considered very accurate as are the precision 
potentiometer and Pailyaneueber used. Since these errors are in excess 
of those guaranteed by the manufacturer, it is recommended that a new 
model be evaluated before any conclusions are drawn as to the 
feasibility of such a component. 

The ability of the sine-cosine potentiometer to null the gravi- 
tational error of a linear accelerometer, as demonstrated in Figure XXX, 
shows promise. The error is consistent with that of the potentiometer 
calibration data and, as was stated previously, a new potentiometer 
should improve the results. 

The integrator response, shown in Figure XXXII, is quite satis- 
factory since an output variation at the recorder of only +2% is ex- 
petienced exer the entire anticipated operating range. The acceler- 
ation signal generated by the Scotch Yoke (Figure X) is itself subject 
to error because of play in its connection to the gear system. This 
could allow slapping at the peak displacements, thus giving an 
acceleration singularity and higher than actual displacement. This 
particular error will not appear in the actual system because the 
mounting will be rigid. 

The frequency response curve shows that the double integrator is 


very good over the entire operating range. The slope is a constant 
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~20 dg per decade which shows the component is operating properly. The 
spread between the theoretical and experimental curves is merely a 
gain that was not predictable in theory and actually works to the 
advantage of the overall system. It is recommended that runs at an 
expanded time scale be made in order that phase differences between 
input and output may be measured. It should be remembered that the 
DeGie amplifier introduces a 180 degree phase shift, and therefore 

the output displacement's phase should correspond exactly to that of 
the input acceleration. 

The expanded Sanborn recording at the beginning of the drift run 
(Figure XXXI) shows the system output in detail. It appears to be an 
undistorted sine wave, and therefore a close reproduction of the input 
displacement to the accelerometer. The noise superimposed on the 
acceleration chart is principally bearing noise from the Scotch Yoke 
and should be reproduced as a slight bumpiness in the displacement 
plot, for the motion is rough. A close examination of the plot will 
disclose this irregularity, thus attesting to the accuracy of the 
system. It should be reemphasized that this type of noise will not be 
present in the actual system. 

The drift of the system is supposed to be negligible for the 
duration of a model run in the towing tank, a period of from two to 
three minutes. The ten minute drift run showed no drift over the 
normal run time and a very slight amount at the end of ten minutes. 
Since the double integrator is an imperfect integrator by intent, it 
should merely amplify and not integrate D.C. and very low frequency 
signals. Therefore, the recorded drift may be due to drift in the 
recording system and not in the integrator. This possibility should 


be investigated further by noting the drift of the recorder at zero 
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input, and then by noting the integrator drift at zero input. 

A further test, suggested by the above discussion, is the response 
of the integrator to a step function. The gain of the integrators at 
zero frequency is approximately 100. Therefore, the output level 
should rise to 100 times the input level and then stop. 

The Sanborn recording of the response of the system to a constant 
displacement at various frequencies shows a random noise in the dis- 
placement plot for 108 rpm. This is the same type of noise ex- 
perienced in runs at a table sme Linea of 45 degrees. The natural 
frequency of the tilt table and its mounted equipment was found to be 
approximately 6 cycles per second. The noise frequency was .16 to .33 
cps at 108 rpm and 0.1 to 0.2 cps at 76 rpm. An investigation should 
be made to determine whether this noise is due to mechanical or 


electrical sources. 
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V. CONCLUSIONS 


The sensing elements selected are believed to be the optimum 
presently available components. 

Space and weight limitations of model carried equipment are 

no longer a problem. 

Model carried components must be positioned and aligned in the 
model with extreme accuracy or serious cross coupling errors 
will result. 

The gravitational error due to rotation of the linear 
accelerometers into, and out of, the earth's gravitational 
field is the only error signal of the system large enough to 
require a special compensation system. 

Servo driven sine-cosine potentiometers are presently considered 
the optimum method for generating gravity compensation signals. 
This is primarily due to their low cost, simplicity, and 
reliability. However, insufficient data has been obtained to 
fully evaluate their effectiveness. 

The proposed system is capable of resolving the low level 
Signals of the accelerometers and performing double integration 
over the entire anticipated frequency range with no noticeable 
attenuation of the output displacement signal. 

The balancing system of the integrating amplifier, at the 
present stage of its development, is too critical for ease of 
operation and should be studied further. 

The drift characteristic of the system is very satisfactory. 
Large amplitude, low frequency, noise present during dynamic 


gravitational compensation studies must be eliminated before 
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further analysis of the linear displacement system can be made. 
In view of the above, the proposed sensing system looks very 
promising. However, more analysis must be done before it can 


be properly evaluated. 


¢ 
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VI. RECOMMENDATIONS 


It is felt that the proposed system is adequate for motion 
measurements but that further studies should be made to refine 
the following components: 

a. The gravitational compensation system. 

b. The integrator balancing system. 

Rate gyros should be procured and used as angular sensing 

elements in order to evaluate their performance while 

measuring angular motions as well as exciting the gravitational 
compensation system. 

Further evaluation of the laboratory mockup should be made through 
the following tests: 

a. Expanding the time scale of the recorder to permit 
measurement of the phase difference between acceleration, 
velocity, and displacement. A phase lag is believed to exist 
Since the amplifier operating range, positioned on the integrating 
slope of Figure XXxXV , is too near the first break point of that 
curve, 

b. Graphical integration of the input signal for comparison 
with the output displacement, in conjunction with variation of the 
Scotch Yoke speed at a forty-five degree angle of tilt, to determine 
the nature of the interference signal superimposed on the dis- 
placement record. 

c. Observation of the response of the Double Integrating 


Summing Amplifier to a step function input. 
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APPENDIX A 


DERIVATION SUMMARY 


1. Expression for the Direct Axis Gravitational Errors: 
a. Heave Accelerometer 
From Figure XXXIV (a): 


* ee x ty = Bg? 


where 

x =e tan ¥ 
and 

y =e tan 9g 


Combining the above-:thfee expressions: 
e@ + eX tan2W + e* tan? g= ge 
e~ (1 + tan?W + tan2 g) 


e~ (sec“¥ + sec* g - 1) 


ge 


2 il 1 2) 
e at 
(55 i cos-g ) : 


ae ce » coseY- cos? VY cos2 )- 2 
=TE 
cos“ Veos“g 





Dividing through by the coefficient of e* and 


extracting the square root: 


ex cos ¥ cos 


ee, 
(cos*¢ + cos~ y _ cos°¢ cos? V’)2 


Eom 
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AUXXX JUNDI4 
A893 OWO19 | dd JBPOWOID | 999 JB JOWOIS | dd 
SIPIS ebans @APOH 
(>) (9) (e) 
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b. Surge Accelerometer 


From Figure XXXIV(b): 


e- + we + 2° = ge 


where 
* y=etan @ 
and 


e cot ¥ 


N 


Combining the above three expressiong and’ °<. 
simplifying in a manner similar to that for the heave 


accelerometer: 


e=g sin ¥cos 6 (2) 
(cos2@ + sin2¥- sin2 ¥cos29)2 


c. Sidle Accelerometer: 


From Figure XXXIV (c): 
ef + x2 4 22 = pA 


where 


e tan 9 


x 
and 


e cot g 


Zz 
Comgining the above three expressions and 
simplifying in the same manner as in part b: 


e=g¢ sing cos 9 


ry 
(cos?0 + sin2% - sin2¢ cose)? (2) 


= Ale 





2. Taylor Series Expansion for the Heave Accelerometer Direct 


Axis Gravitational Error: 


ess cos ¥ cos g 
(cos? + cos ¥ -cos* g cos“ ¥) Cy.) 


cos Ycos g = 3 | cos( Ws J) + cos(¥- g) | 


Ss 


(cos coe g)- = 2 | cos*( w+ 3) + cos2(¥- 9) + 
2 cos( + Z) cos( ¥- g)] 


cos? ¥ cos2g = 


Pe 


| cos2( Ys BS) + cos<( Y/. B)+ cos 2 scos 2g] 


cos? ¥ = i (1 + cos 24) 
cos*G = (1 + cos 2 Q) 


cos?( by G) = 4/1 + cos 2 (%+ ) | 
cos“(W- g) = ah + cos 2 (¥- 2) | 


Substituting the above relations into equation (1): 


one. ____ Eleosty+m+costy-pJ 
g ($+4c0s Lfr+zrgcos2y mar -g cos 2(y+g)-¢ —L-cos &(y-9)- Fcos2y-tcos 2p) 


+z[cos2p+cos 2y ~Scos2ly- p)-Feos Bip 


@ 
! 
0a 


a _ cos(y+p)+cos(y-f) 


pi cos CMa (Ja) 
* [3+cos ed +C0S2y -5C0S 2(¥-f)-sc0S a(p+ | * 


Eo 





Expanding Eq. (1a) into a Taylor Series: 

: V+ 2yd +g  » yrs ay? 62 Ades bt 

cos( + J) =I- : + Ve 4d +6 YO 4Ve + or 
Z4 
"6b +15 444 20747 4 sytdtequdes df 
720 
cos(f =o) =|— a Ee 1 = Th +6u'o*- AY oat 4 
_ 8-645 + 15 44¢7- 2047674 15 gto g74 oF 


120 


Summing the above gives an expression for the numerator of (1a): 


4 
Num. = 2-Y—$ + we ei ue - Hee Aree eo 


360 
Expanding the denominator terms of (la): 
ie Z 2 44 6 6 
Cos 26=| 2p += 4 KpP+e-- 
es te 14 eee ee 
Cos2h= 1-245 +— yp -S y+ 


-FCos2ly-N—=-5 4-249 4 97-4 (Vt ayros o¥Oe 409% 04) 
+& (V6 y+ 15 Y'g-20y'd'+ 5Y'O*-6YG°4 O°) +- = - 
-1 Cos 2(4+4) =-5 ++ 24b4 0-3 (44 474 6d + 440? + O° 
+5 (yo+ 6 y'$ +15 4074 2o¥ G+ 1529's 64074 go) - * 
Combining terms the denominator becomes: 
I 
Denom. = (4- 4y'q? + : [ytgt+ ye gt]y.--)2 


The series expression for ( / ) Ls, therefore: 


og tie P te VO + (YEO) ra (U"4 V2 sen (WO4. GP) pee | 
: (4-447 74 5 (vig sytpt]+.- -y2 (1b) 
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3. Derivation of the Operational Amplifier Feedback Loop 

A schematic diagram of an electronic integrator is shown in 
Figure XXXV. A perfect integrator would not have Rp across C. 
Re was placed there to prevent the integration of D.C. and low 
frequency signals, thus reducing drift. It is hoped that the low 
gain at high frequencies will prevent interference from high fre- 


% 
quency noise, 


Oo. = alt 
i -A or Cr = x %o (1) 
g 
e; - e, = iRy (5) 
e,-¢e = -e,(441) = 12 
g o> ORK f 
where 
Re 
12> 1.10-8 
7 10 
Therefore 
Ce ee a 
l + jwCRe 
i = -e,(1 Se jwRe) Ca) 
Re 


Substituting Eq. (4) and (6) into (5) 


0 
es +B - “Eo Ri (Gil + jooCRe) 
A Re 
Therefore 
Go _ = R 
el : (7) 


a 
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FIGURE XXXV 








Plot (a) of Figure XXXV makes use of a real time integrator 
(R;C = 1 sec) modified by Rp. The frequency to be used in the mockup 
(Figure X) is approximately 8 radians per second (76 RPM). At this 
frequency the gain is ~9.15 dg or 0.122. This is not very satisfactory 
for a gain of at least unity (for recording purposes) should be 
attempted. 

It can be seen from Figure XXXV that if the plot were to cross 
the O dg. line at a frequency of 10 radians per second, the time 
constant (CR,) would have to be 1 second for a steady state gain of 10. 

If C were to remain 1.0 microfarad , Rp would have to be 
reduced to one megohm. It is desired to keep R; at least one megohm., 
Therefore, the steady state gain would be reduced to 1 (0 dg). This 
is illustrated in plot b. This would gain nothing. 

However, if the capacitance is varied, the steady state gain 
remains constant, and only the break point varies. Let the steady 
state gain be 10 (Rp = 10 R; = 10 megohms). To have gains greater 
than unity over the operating range, let C be 0.1 microfarads. The 
gain at calibration frequency will then be +1 dg or 1.25. This is 


satisfactory. (See plot (c)). 
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Colum #1 
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20 
30 
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Column #1: 
Column #2: 
Column #3: 


Colum #1: 


APPENDIX B 


Summary of Data and Calculations 


TABLE III 


Calibration of Accelerometer No. 6429 


s #2 


32.105 
21076 
BO) acl 
27.855 
21.638 
20.676 
16.083 
11.000 
558k 
0.000 

- 5.56), 
-11.000 
-16.083 
20.676 
—21,.638 
27.855 
-30.225 
-31.676 


#3 


17.0 
16.7 
15.93 
14.671 
ey, 
10.854 
8.39h 
5 WLS 
2.868 


= 40525 


- 2.978 
- 5,825 
- 8.480 
-10.880 
-12.938 
-14,.612 
~15.831 
-16.58 


#l, 


3.64 

2eGH 

2.49 

Zl 

Dey) 

0.94 

-0,21 
-0.65 
-0.98 
71.00 
1s 

1.45 

1.44 

1.44 

ale 1s: 

1.01 

sy 

- .57 


Angle of Rotation 


Acceleration in fps 


2 


Output in Millivolts 


Column 


180 
190 
200 
210 
220 
230 
210 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 


Non-linearity and hysteresis 


“772 


#1 #2 


—32.165 
-31.676 
—30.225 
27.555 
-21,.638 
—20.676 
-16.083 
-11.000 
>a Pager 
0.000 
558k 
11.000 
16.083 
20.676 
24,.638 
PH Lolo 15) 
B0ee25 
Bio/0 
32.165 


error as a % of z(max) 


#3 


-16.80 
-16.54 
-15.755 
-14,.498 
-12.802 
-10.736 
- 8,326 
- 5.670 
- 2.82 
+ 0.066 
2.985 
5.808 
8.457 
10.847 
12.899 
14.565 
15276 
Gm: 
16575 


th, 


- 19 
- 19 
-1.05 
= 
-—1.44 
-1.32 
1.51 
-1.51 
-1.47 
71.26 
1.26 
leg 
1.00 
0.80 
0.42 
Sha eks: 
-0,.36 
-0.77 
-1.15 





TABLE IV 


Calibration of Accelerometer No. 3023 


Column #1 = #2 #3 #1, Column #1 #2 #3 #1, 
O 0.000 0.080 -8.27 180 0.000 02535 + 6.37 
10 § 584 2.750 -8.62 190 — 5epeh = 25159 = 5295 
20 11,000 5.396 -8.97 200 -11.000 - 4.779 - 5.28 
30 16.083 72801 -8.20 210 -~16.083 = 7235 —- 4.73 
L0 20.676 10.033 ~7.30 220 ~20.676 = 92163 = 3.96 
50 214.638 11.974 -6.09 230 -21,.638 -11.389 —- 3.22 
60 2fmeo> “182560 =h.73 2h0 —27.855 -12.957 - 2.44 
70 30.225 Lig ya —3.41 250 =30.225 -liei2 = isc 
80 31.676 15.489 -1.64 260 -31.676 -14.848 - 0.58 
90 32.165 15.775 0 270 —32.165 -15.102 O=52 
100 31.676 15.585 1355 280 —-31.676 -14.902 NAL 
110 S0g225 14.941 3.92 290 ~30.225 -—14.239 Zeee 
120 27 855 13.840 L428 300 —27.855 -13.141 3.48 
130 24.638 12.337 DO 310 -21,.638 -11.643 1.96 
140 20.676 10.462 6.50 320 -20.676 - 9.796 6475 
150 16.083 8.276 7.08 330 -16.083 —- 7.633 8.08 
160 11,000 5.819 6.50 340 -11,.000 —- 5.239 955 
170 5 58h 3.226 6.70 350 ~ 5.584 - 2.668 1Om41L 
360 - 0,000 - 0.004 10.96 
Colum #1: Angle of Rotation 
Column #2: Acceleration in fps* 
Column #3: Output in millivolts 
Column #4: Non-linearity and hysteresis error as a % of y(max) 
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Column #1 


0 
10 
20 
30 
L0 
50 
60 
70 
80 
90 

100 
110 
120 
130 
140 
150 
160 
70 


Colum #1: 
Column #2: 
Column #3: 


Colum #4: 


TABLE V 


#2 #3 #l, 
0.000 - 0.041 -6.79 
5 58h o> 9 —O.05 

11 600 5.441 -6.70 
16.083 8.015 -6.02 
20.676 10.346 -5.32 
24.638 12.389 ~4.25 
27.855 14.045 -3.40 
BOeee5 15.290 =-2,38 
salyotde 16.074 -1.43 
B2nroS 16.41 0 
31.676 16, 2k IL AL 
30.225 15.578 2.06 
212099 Tish7O 3616 
21,.638 12.941 4.27 
20.676 DAL AOUES: 1.96 
16.083 Seyi Sis Sy 
11.000 6e2)1. 75219 
5 58h ols wees 

Angle of Rotation 

Acceleration in fps 

Output in Millivolts 


Non-linearity and hysteresis error as a % of x(max) 


Calibration of Accelerometer 


Column #1 


180 
190 
200 
210 
220 
230 
2h0 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 


~19= 


No. 3024 


#2 


0.000 
= Sou 
-11.000 
-16.083 
-20.676 
—21,.638 
—27.855 
-30.225 
-31.676 
—32.165 
-31.676 
—30.225 
27.855 
-21,.638 
-20.676 
-16.083 
-11.000 
a 5.584 

0.000 


#3 


0.703 
Seren eal 
- 4.824 
- 7.371 
- 9.696 
-11.696 
-13 .330 
-14.547 
-15.309 
-15.612 
15.427 
-14.782 
-13.689 
—25177 
-10.285 
- 8.071 
= ae 
- 2.985 
- 0.164 
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TABLE VI 
Evaluation of Exact Expression for Gravitational Error 


in the Heave Accelerometer at Various Angles 
of Roll for Fixed Angles of Pitch 


(°/g) exact; Eq. (1) 
Y=? Y= 5° Y= 10° 


s 0 1.000 9962 98L8 
10 9848 9812 ~9'703 
20 Sei 9366 29270 
30 8660 8636 8559 
L0 «7660 7643 «1589 
50 6428 6418 6387 
60 . 5000 AINSI 4981 
TABLE VII 


Percent Difference Between Gravitational Signal of the 
Heave Accelerometer for Various Pitch Angles 
Expressed as a Percent of the Maximum 
Predicted Heave Signal (10 fps) 


“% Diff. between 
columns of constant W 


Go 
0° - 5° 0° ~109 5° ~ 100 

0 1.22 4.89 3.61 
10 1516 Le67 Saou 
20 1,06 4.09 3.09 
30 Te 34> Zee 
1,0 aoe 2.29 Pewee 
50 032 1432 1.00 
60 16 <6 045 


¥ % Difference is referred to the maximum heave signal 
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go 


10 
= 
20 
25 
30 
BD 
L0 
L5 
50 
DD 
60 


TABLE VIII 


Percent Error due to Series Approximations of 
Gravitational Error in the Heave Accelerometer for 
Various Roll Angles at Constant Pitch 
Angles of Zero and Ten Degrees 


fe in percent of maximum heave signal 


. Y=0 ¥ = 10° 
. (1b) Eq. (1b)* Eq. (1b)™" Bq. (ab) Eq. (1b)* Eq. (1b) 

O O 0 O 0 0 

0 0 0 0322 0 2576 
O O 0 0966 O 02576 
0 0 0322 O - 06h, STL 
0 .0322 1932 0322 .0966 1.030 
O 0 483 .1610 0 1.835 
0 O 9982 -.0644 0644 2.673 
0 O 1.868 - .0966 - .4508 3.928 
0 0644 3.091 =0222 = 45 h2 5.732 
O .0966 hago +.0322 - .4830 8.179 
0 1932 7.599 Mizss’ = sea72. | lea 

O 03542 11.045 6225) 48=—- = eee 14.97 

0 ©5796 553 sone ~1.739 19.83 


* Eq. (1b) through fourth order terms. Appendix A. 


** Eq. (1b) through second order terms. Appendix A. 


ae 2 i = 





TABLE IX 


Analysis of Gravitational Error in Sidle and Surge Accelerometers 


g #1 #2 #3 wy # #5 


5 0872 .0871 
10 eer e736 
15 62588 =. 2587 


.0537 5 .0871 -.0250 
0537 10 .1736 =5@250 
0537 15 4.2587 see250 


20 .3420 .3419 = .0537 
25 We ae een223) 06 = 1610 
30 -5000 .4995 — .2684 
20 -Difom 5728 = 2120p 
1,0 6428 6418 = 5367 
45 S707 16 (@5/ = = a (5K 
50 .7660 .7642 - .9660 


> See oly =) 7127 
60 .8660 .8635 -1.342 


Column #1 and #2 are Eq. (3) Appendix A divided by "g" 
evaluated at 6 = 0° and 5° respectively. 


Column #3 is the percent difference between columns #1 and #2 
referred to a maximum sidle signal of 6 fps@. 


Column #1 and #4 are Eq. (2) Appendix A divided by "g" 
evaluated at 9 = 0° and 5° respectively. 


Column #5 is the percent difference between columns #1 and #4, 
referred to a maximum surge signal of 12.5 fps¢. 
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TABLE X 


Calibration of Gamewell RL-~-L4MS SineCosine Potentiometer 


Colum #1 


OO O~AIRNWN EWN FF 


10 
val 
IW 
a5 
14, 
aby 
16 
19 
18 


Column fl: 
Column #2: 


Column #3: 


#2 


OisO5 
MOLL es: 
Dive6> 
Wz ie 
USL GOS 
141.65 
151.65 
161.65 
ESOS 
ielo> 
191.65 
201.65 
ZA Ale: 
221.65 
25 (oes 
ZANE 
251.65 
261.65 


#3 


63 
38 

- oll 
098 
69 
eee 
deg 
126 
1.99 
16 

- ff 
IL ole: 
64 
B85 
095 
02 
04 

- .69 


Column #1 


Rotational Position 


Pot. Setting in Degrees 


332 


19 
20 
Zaih 
22 
23 
2h 
25 
26 
a7 
28 
29 
30 
Sal 
32 
Be, 
34 
aD 
36 


Error as a % of maximum signal 


#2 


Zyl els 
281.65 
291.65 
301.65 
Bile ou 
3212 Gr 
BSL so 
341.65 
Blaser 
1.65 
11.565 
Zales 
31.65 
41.65 
planes 
61.65 
71605 
ells 





Column #1 


Column #1: 
Column #2: 
Column #3: 
Column #1: 


Column #5: 


TABLE XI 


Ability of Sine-Cosine Potentiometer to Null 


Se 
1,0 
% 
ce 
30 
25 
20 
15 
10 


y) 
0 


Gravitational Errors 


#2 


15.236 
13.890 
12.180 
iC) eal 
9.205 
7 «480 
5.655 
3.735 
1.908 
0.000 


#3 


Leeen6 
13.850 
12.360 
10.664 
8.964 
12250 
5 416 
3430 
1.606 
- .020 


Accelerometer Output (mv) 
Potentiometer Output (mv) 
Column #2 minus #3 


Percent Error Referred to "g" 


~8),— 


th, #5 
O 0 
04 ~186 
a2 0097 
149 opt 
2h1 Pate 
Aaah, 07 
e207 I Tal 
0305 1.42 
e902 1.40 
.02 093 


Inclination of Accelerometer from Horizontal 





TABLE XII 


Frequency Response of Double Integrating Amplifier 


RPM 108 76 56 hs} 15 
Rad/sec ia, 3 1298 5.86 3.46 Deas (i 
Cyc¥es/sec. 1.80 1.27 933 5 025 
Acc. Input (mv)” 16.65 8.39 4.53 1.58 325 
Output (mv) 2950 3000 2950 2800 2125 

G= ——_. 1.772 3.58 Ges ve eigae nos. 
dg = 10 log = |G| 2.50 5.54 8.04 12.50 £18.17 


F = Calibration Factor (Accelerometer No. 3024): 2.571 mv/volteg 


* .', acceleration input = Vy x a 


where a = 23 tt = Seg x 47 2(RPM)* 
: a 


a = .0213 (RPM)@ 


A = Video amplifier gain = 100 


Vp = Accelerometer excitation voltage = 8.04 volts (DC) 


Zor 





APPENDIX C 
Sample Calculations 


1. Calculation of Tangential Accelerations Due to Roll, Pitch, 


and Yaw 


a, = (axial separation ee or 


“Op = Ue ins) (2,88 rad/sec~) = .96 fps 
aby > Bind) (1.92 rad/sec*) = 1.28 fps® 
at n = (Sanu) (.54 rad/sec2) = .36 fps 


2. Calculation of Linear Accelerometer Errors Due to Rotational 


Metions @f themModel.” 







yer 


fF, 5 ny BP 2.52 7 | ay 2.52 r 
= - - = 
Te 2 eS 


“y ~ “¥g * “Yq 


ae = ele = 
r 4 a Ss o 6r 
ay = 13% — O45 le y 12 Y Bll 


* See Figure XXXVI 
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Calculation of Resistance Values to be Used in Sine-—Cosine 


Potentiometer Voltage Divider Network 





SLne - Cosine 
1 Potentiometer | 


1) E, varies sinusoidally from O—>-|E,| 


2) Consider adding Ry, 


BE. -E, -E I EL 
p= = a S=o 
? Roy & Roo Ry 
Bee UE 


Ey Rp, (Rpo + Ry) 


Ry (Rpg + RL) + Ryy By 


= mi Rpg Ry, | 


Let R, >» Ry, 3 Rr) Rp, 


ae En = E | Apo =| 
ss tia Roo 


~88— 





4, Calculation of Motor Speed to Obtain a Given Acceleration, ata 


Particular Displacement Amplitude 


a | s | Sinw t |s| - 4.76 — ~ litude) inches 
G) = = rad/sec 
Differentiating 
a= = x Lf Sinwt Let [al = 12.5 fps? 
af 


2 
Pe pe = -4.763 (Lif) me = 628 


IEA T= .792 sec/cycle 


f = 1_ = 1,262 cps 
792 i 


RPM = 60 (1.262) = 76 


5. Calculation of Signal Attenuation Upon Integration 


Consider a=A Sinw&wt A is acceleration amplitude. 


v=z> * CosWt 


o 
i 


2 
- Gz Sinwt or s=-A[ 7 SinWD t 


if = g lL: Signal voltage will be attenuated upon integration 


J => 1: Signal voltage will be amplified upon integration 
Therefore, in order for the signal to be amplified, the period must be 
greater than 6.28 sec. Since the predicted periods are approximately 


one second, the integrating circuits must have an amplification of at 


least 6.28 to prevent signal attenuation. 


accu 





APPENDIX D 


Equipment Used in Various Tests 


1. Static Calibration of the Linear Accelerometers 
(a) Hardinge Milling Machine Dividing Head, Serial #705. 
Hardinge Brothers, Inc., Elmira, N. Y. 
(b) Four jaw lathe gee, five-inch, Mo. 9405-S, Serial #44. 
Union Mfg. Co., New Britain, Conn. 
(c) Rubicon Precision Potentiometer, Serial #42765. 
(d) Rubicon Galvanometer (external) for use with item (c), 
Serial #E-062. 
(e) Vacuum Tube Voltmeter, Simpson Electric Company, Chicago, 
Model 303. 
(f) One Burgess F4BP 6.0 volt battery. 
2. Static Calibration of the Sine-Cosine Potentiometer 
(a) Items (a) through (f) of paragraph number 1. 
(b) Sine-Cosine Potentiometer, Gamewell RL-14MS, Serial #674, 
Terminals B-D, Card Resistance 35 , L00 ohms. 
(c) Two "Ohmite"” resistors, 3.6 meg-ohm, and 10,000 ohms 
respectively. 
3. Ability of the Sine-Cosine Potentiometer to Null Gravitational Errors 
(a) Items (c), (a), and (e) of paragraph l. 
(b) The calibrating table shown in Figure IV. 
(c) Statham A6a~1.5-350 linear accelerometer Serial #3023. 
(d) One Burgess 4FH 1.5 volt battery. Potentiometer excitation. 
(e) One Burgess 5540 7.5 volt "c" battery. Potentiometer excitation. 
(f) Decade Resistance Type 602-J, Serial #5041, General Radio, 


Cambridge, Mass. 
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(zg) One 3.6 meg—ohm "Ohmite" carbon resistor. 


(h) 
(i) 


Item (b) of paragraph 2. 
One Willard DD-~3-3, 6 volt, three cell, automobile type 


battery. Accelerometer excitation. 


(j) One Willard DD-7-1, 2 volt, single wet cell, battery. 


Accelerometer excitation. 
s 


Double Integrator Response to a Constant Displacement Input 


(a) 


(b) 


(c) 


(d) 


(e) 


Philbrick Power Supply with Sorensen Regulated A.C. Power 
Supply. Sorensen & Co., Inc., Serial 2+664, Model 2000-S. 
The Double Integrating Summing Amplifier shown pictorially 
in Figure III, and schematically in Figure V. 

Twin-Viso Recorder, Sanborn Company, Cambridge, Mass. 
Model 60-1300, Serial #24. D.C. Pre-Amplifier Model 60-400, 
Serial #114. D.C. Amplifier channel "A" Model 64-300A, 
Serial #393, channel "B" Model 64-300A Serial #376. 

D.C. Amplifier, Video Instruments Co., Inc. 2340 Sawtelle 
Blvd., Los Angeles 64, California. Model 71, Serial 129. 
scotch Yoke Drive, see Figure X. 

Motor Data: 


Electro Machines Inc., Cedarburg, Wisconsin. 


A.C. Motor Volts 220 
Type P Amps L 
Frame D42S 60 cps 

H.P. 1.0 3 phase 
RPM 31,50 Rating 40°C 


Variable Speed Transmission 


Graham Transmission Inc., Milwaukee, Wisconsin. 
Model 225 
Serial #73039 
Output RPM 1200 
Input RPM 3600 


(f) Items (b), (c), (i), and (j) of paragraph 3. 


(g) Item (e) of paragraph 1. 


Bone 





5. Dynamic Gravitational Error Compensation of the Linear Accelerometer 
Signal at Fixed Angles of Inclination 


(a) All of the equipment listed under paragraph 4. 
(b) Item (b) of paragraph 2. 
(c) Items (d) and (e) of paragraph 3. 


(d) Decade Resistor 11.1-k ohm, General Radio. 
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